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Abstract 

This paper presents updated measurements of the lifetimes of the Bg meson and 
the AjJ baryon using 4.4 million hadronic Z° decays recorded by the OPAL detector at 
LEP from 1990 to 1995. A sample of Bg decays is obtained using D~£ + combinations, 
where the D~ is fully reconstructed in the 4>tt~ , K*°K~ and K~Kg decay channels and 
partially reconstructed in the 4>£~uX decay mode. A sample of A° decays is obtained 
using A£l~ combinations, where the A+ is fully reconstructed in its decay to a pK~7r + 
final state and partially reconstructed in the A£ + zaK decay channel. From 172 ± 28 
D~£ + combinations attributed to Bg decays, the measured lifetime is 

r(B s ) = 1.50l°i|±0.04ps, 

where the errors are statistical and systematic, respectively. From the 129 =L 25 Af£~ 
combinations attributed to A^ decays, the measured lifetime is 

r(A°) = 1.29t8jf±0.06ps, 
where the errors are statistical and systematic, respectively. 
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1 Introduction 



The lifetimes of b flavoured hadrons are related both to the strengths of the b quark couplings to c 
and u quarks, described by the CKM matrix elements V c b and V u h, respectively, and to the dynamics 
of b hadron decays. The spectator model assumes that the light quarks in b and c hadrons do not 
affect the decay of the heavy quark, and thus predicts the lifetimes of all b hadrons to be equal. 
For charm hadrons this prediction is inaccurate; non-spectator effects, such as interference between 
different decay modes, result in a D + lifetime approximately 2.5 times that of the D° and more 
than twice that of the D~ Q. Models which attempt to account for non-spectator effects predict 
that the differences among b hadron lifetimes are much smaller than those in the charm system due 
to the larger mass of the b quark ]2|, ||, |4|. These models predict a difference in lifetime between 
the B + and B° meson of several percent, and between the Bg and B° meson of about 1% g 1- 
OPAL ||, |j, and other collaborations P|, have published measurements of the Bg lifetime which 
are in agreement with these models. 

Non-spectator decays of B mesons proceeding via W-exchange are Cabibbo-allowed but are ex- 
pected to be suppressed, relative to spectator decays, by an amount depending on the ratio of 
the initial-state meson mass and the final-state quark masses (helicity-suppression). This suppres- 
sion does not occur for baryon decay, therefore b baryon lifetimes are expected to be shorter than 
b meson lifetimes. This expectation is consistent with existing lifetime measurements [pj for which 
t(A^)/t(B°) = 0.73 ± 0.06 Q. It is also supported by a recent theoretical prediction ||] which 
yields r(Ag)/r(B°) of about 0.9. Reference f§ finds this ratio to be 0.98 + 0(l/m|). Measure- 
ments of the average b baryon lifetime have been published based both on analyses of A£~ and 
A^£~ correlations by OPAL |^, [| and other collaborations [JO]]. In both analyses, the dominant 
contribution is expected to come from A{J baryons, though both Al~ and A££~ combinations can 
arise from the decays of other b baryons. The composition of each sample depends on the b baryon 
production fractions, but the A^£~ correlations provide a purer sample of A^ baryons. 

This paper presents updated measurements of the Bg meson and A° baryon lifetimes using D~^ + 
and A^£~ combinations reconstructed from the full OPAL hadronic data sample collected on or 
near the Z° resonance. These results supersede the previous OPAL measurements using D~^ + || 
and A~^£~ Q combinations. The decay channels used for these lifetime measurements are:[| 



where i is an electron or a muon. In each case, the proper decay time of the b hadron is determined 
on an event-by-event basis using measured decay lengths and estimates of the b hadron energy. 

The following section provides a brief description of the OPAL detector. The remaining sections 
describe the selection of D~£ + and A^£~ candidates, the determination of the b hadron decay 
lengths, the estimation of the b hadron boost, the lifetime fits, the results, and the systematic 
errors. 



^dn this paper, charge conjugate modes are always implied. 



B° -» D~ i + v X 

U K*°K- 




K*° -> K+vr 
(j) -» K+K 

Kg -» 7T+7T- 

<p -> K+K~ 



A° -> A+ e~vX 

' — ► pK _ 7T + 

U M+vX, 



A 



P7T 



U (f)£-i?X, 



1 



2 The OPAL Detector 



The OPAL detector is described in detail in reference The central tracking system is com- 

posed of a precision vertex drift chamber, a large volume jet chamber surrounded by a set of 
chambers which measure the z-coordinateQ and, for the majority of the data used in this analy- 
sis, a high-precision silicon microvertex detector. These detectors are located inside a solenoidal 
coil. The detectors outside the solenoid consist of a time-of-flight scintillator array and a lead 
glass electromagnetic calorimeter with a presampler, followed by a hadron calorimeter consisting 
of the instrumented return yoke of the magnet, and several layers of muon chambers. Charged 
particles are identified by their specific energy loss per unit length, dE/dx, in the jet chamber. 
Further information on the performance of the tracking and dE/dx measurements can be found in 
reference [12|. 



3 Monte Carlo Simulation 



Monte Carlo simulation samples of inclusive hadronic Z° decays and of the specific decay modes of 
interest are used to check the selection procedure and lifetime fit procedure. These samples were 
produced using the JETSET 7.4 parton shower Monte Carlo generator [13| with the fragmentation 
function of Peterson et al. [14] for heavy quarks, and then passed through the full OPAL detec- 
tor simulation package [p^|. A special sample of simulated data was generated using a modified 



JETSET decay routine for b baryons [16], where it is assumed that the polarization of the b quark 
is carried by the b baryon. An additional form factor [jlT]] describing the energy transfer from the 
b to c flavoured baryon was used in the generation of the polarized sample. 



4 Candidate Selection 



This analysis uses data collected during the 1990-1995 LEP running periods at centre-of-mass 
energies within ±3GeV of the Z° resonance. After the standard hadronic event selection [18] and 



detector performance requirements, a sample of 4.4 million events is selected. Jets are defined 
using charged tracks and electromagnetic clusters not associated with a charged track. These are 
combined into jets using the scaled invariant mass algorithm with the E0 recombination scheme [|l9| 
using y cut = 0.04. 

Only charged tracks that are well-measured in the x-y plane are considered in this analysis. Well 
measured tracks are defined according to standard track quality cuts pCfl . Within a single jet, not 
all combinations of accepted tracks with the appropriate charge combination are considered in the 
D;T and A+ searches. Instead, to reduce the combinatorial background, the dE/dx probability, 
wi, that the observed dE/dx is consistent with the assumed particle hypothesis, i, is required to 
be greater than 1%. More restrictive requirements are imposed on a channel-by-channel basis, as 
described in the following subsections. 

Leptons are identified as follows. Electron candidates with a momentum of at least 2GeV/c are 
identified using an artificial neural network based on twelve measured quantities from the elec- 



2 The right-handed coordinate system is defined such that the z-axis follows the electron beam direction 
and the x-y plane is perpendicular to it with the a:-axis lying approximately horizontally. The polar angle 9 is 
defined relative to the +z-axis, and the azimuthal angle cf> is defined relative to the +x-axis. 
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tromagnetic calorimeter and the central tracking detector pl| . Between 1 and 2 GeV/c, electron 
candidates are required to have a dE/dx probability of larger than 1% for the electron hypothe- 
sis and less than 1% for the proton hypothesis, because it is in this region that the electron and 
proton dE/dx bands cross. Electron candidates identified as arising from photon conversions are 
rejected [E2]. Muons are identified by associating central detector tracks with track segments in 



the muon detectors and requiring a position match in two orthogonal coordinates [22 



4.1 Selection of D s candidates 

The D~ candidates are reconstructed in four modes: 

1. D~ -> K*°K~ in which the K*° decays into a K+vr". 

2. DjT — > (j)TT~ where the (j) subsequently decays into K + K~. 

3. DjT — > K~Kg where the Kg decays into it + it~. 

4. The D~ is partially reconstructed in D~ — > c/rf~~PX where the <fi decays into K + K~. 

In all modes except (J)£~vX, charged kaon candidates are required to have a momentum greater than 
2GeV/c. Also, because of the potential for misidentifying a pion as a kaon, if the observed energy 
loss of a kaon candidate is greater than the mean expected for a kaon, it must satisfy wk > 5%. 

In both the D~ — > K*°K~ and the D~ — > (j)Tr~ channels, if the observed energy losses of both 
kaon candidates are greater than the mean expected for a kaon, the product of the two dE/dx 
probabilities must satisfy wki • wk2 > 0.02. The momentum of the K + K~7r~ combination is 
required to be greater than 9GeV/c for both channels. 

For the K*°K~ mode, the invariant mass of the K + 7r~ combination is required to satisfy 0.845 < 
n^Kir < 0.945 GeV/c 2 . To reduce the possibility of mistaking a D~ — ► K*°ir~ decay for the desired 
signal, the measured dE/dx of the K~ candidate must be at least one standard deviation below 
the mean dE/dx that is expected for a pion. 

In the (p7r~ mode, the observed 4> width is dominated by detector resolution and the K + K~ mass 
is required to satisfy 1.005 < itikk < 1.035 GeV/c 2 . The momentum of the (j) candidate is required 
to be greater than 4.0 GeV/c. 

Differences between the angular distributions of D~ decays and those of random track combinations 
are used to suppress further the combinatorial background. The D~ is a spin-0 meson and the final 
states of both decay modes consist of a spin-1 (</> or K*°) meson and a spin-0 (tt~ or K~) meson. 
The D~ signal is expected to be uniform in cos# p , where 9 P is the angle in the rest frame of the 
D;T between the spin-0 meson direction and the D~ direction in the laboratory frame. However, 
the cos 9p distribution of random combinations peaks in the forward and backward directions. It is 
therefore required that | cos9 p \ < 0.8 (0.9) for the K*°K~ (4>tt~) mode. The distribution of cos^, 
where 6 V is the angle in the rest frame of the spin-1 meson between the direction of the final state 
kaon from the decay of the spin-1 meson and the D~ direction, is proportional to cos 2 9 V for D~ 
decays. The cosB v distribution of the random K^K it combinations in the data is, however, 
approximately flat. Therefore it is required that | cos#„| > 0.4. 

For reconstruction of the D~ — > K~Kg channel, accepted charged kaons are combined with Kg 
mesons reconstructed in their decays to ir + ir~ as described in The ir + ir~ mass is required to 
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satisfy 0.475 < < 0.525 GeV/c 2 . The background of Kg particles from fragmentation is reduced 
by requiring the Kg momentum to be greater than 3GeV/c. To reject background from A decays 
where the proton is misidentified as a pion, K g candidates are rejected if 1.10 < m p7r < 1.13 GeV/c 2 , 
where m p7T is the invariant mass of the two tracks when the highest momentum track is assigned 
the proton mass. To improve further the mass resolution of the D~, a fit using kinematic and 
geometrical constraints is performed. The mass of the two tracks forming the K g is constrained to 
the known K s mass [Q]. Further constraints are applied to the D~ and Kg, in which the directions 
of the vectors between their production and decay points are constrained to be the same as the 
reconstructed momentum vectors. 

For the D~ — > 0£~PX selection, the purity of the kaons from the <j) decay is enhanced by applying 
additional pion rejection. This requires that the dE/dx probability for the pion hypothesis be less 
than 40% for tracks whose momenta are such that the mean dE/dx for pions is higher than that of 
kaons, and less than 1% for low momentum tracks where the mean dE/dx for kaons is higher than 
that of pions. Additionally, the product of the dE/dx probabilities for a pion hypothesis for the 
two kaon candidates must be less than 0.018. The <j> candidate momentum is required to be greater 
than 4.6GeV/c. The lepton candidate is required to have a momentum greater than 1 GeV/c for 
electrons and 2GeV/c for muons. The invariant mass of the (f>£~ combination must be less than 
1.9GeV/c 2 . 

4.2 Selection of A+ candidates 

The A+ candidates are reconstructed in two modes; A^~ — > pK~7r + and — > A£ + ^X. In the 
pK~7r + selection, the proton, kaon and pion momenta are required to be greater than 3, 2 and 
1 GeV/c, respectively. If the observed energy loss of the proton or kaon candidate is greater 
than the mean dE/dx expected for that particle, the corresponding dE/dx probability is required 
to be greater than 3%. The dE/dx probability of the pion hypothesis for the proton candidate 
is required to be less than 1%, which substantially reduces the combinatorial background. The 
pK~7r + combination must have a momentum greater than 9GeV/c. To reduce the possibility of 
mistaking the decay D+ — > i^(K + K~)7r + for a A^ decay by misidentifying the K + as a proton, 
candidates are rejected if the invariant mass of the pK~ combination, when the proton candidate 
is assigned the kaon mass, is within 10MeV/c 2 of the nominal <fi mass [jy]. 

In the A+ — > A£ + i/X selection, A candidates are identified via the decay A — > p-7r - . The selection 
procedure is similar to the one used in reference ||. The track with the larger momentum is 
assumed to be the proton and its momentum is required to be larger than 3.0 GeV/c. The other 
track is required to have a momentum larger than 0.8GeV/c. The selection criteria for the lepton 
are the same as described in the previous section for the selection of D~ semileptonic decays. The 
invariant mass of the A£ + combination is required to be less than 2.2GeV/c 2 . 

4.3 D<T^ + ano - ^-t^~ selection and decay length determination 

Once a combination of tracks that satisfies the D~ or A+ candidate selection is found, a search 
is performed to find a lepton from b hadron decay of opposite charge in the same jet. Lepton 
candidates are identified as described in the introduction to section Both the electron and muon 
candidates are required to have a momentum greater than 2GeV/c except in the D^T — » K _ Kg 
mode where a higher momentum cut of 5 GeV /c is used to improve the signal to background 
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ratio. It is also required that the lepton candidate track be measured precisely by either the silicon 
microvertex detector or the vertex drift chamber. 

To further suppress combinatorial background, requirements are made on the invariant mass and 
momentum of the ~D~£ + and A^£~ candidate combinations. K + K~ir~£ + combinations are required 
to have a mass between 3.2 and 5.5 GeV/c 2 and momentum larger than 17GeV/c. K _ Kg£ + combi- 
nations are accepted if they have an invariant mass between 3.4 and 5.5GeV/c 2 and a momentum 
greater than 17 GeV/c. The <p£~£ + mass is required to be less than 4.8 GeV/c 2 and its momentum 
greater than 12GeV/c. Also, for the <p£~£ + mode, the invariant mass of the cj> and the £ + must 
be greater than 2.1GeV/c 2 . In conjunction with the invariant mass cut on the cj)£~ pair, this un- 
ambiguously separates the leptons from the Bg and D~ decays. The pK~7r + ^~ combination must 
have a mass between 3.5 and 5.5GeV/c 2 and momentum greater than 17GeV/c. In the A£ + £~ 
channel, combinations are accepted if the A£ + £~ invariant mass is greater than 2.5 GeV/c 2 and less 
than 5GeV/c 2 and the invariant mass of the A and £~ is greater than 2.2GeV/c 2 . Furthermore, 
the cosine of the opening angle between the lepton and the D~ or A+ candidate must be greater 
than 0.4. 

Three vertices — the beam spot, the Bg (A^) decay vertex and the D~ (A+) decay vertex — are 
reconstructed in the x-y plane. The beam spot is measured using charged tracks with a technique 
that follows any significant shifts in the beam spot position during a LEP fill pif| . The intrinsic 
width of the beam spot in the y direction is taken to be 8 /im. The width in the x direction is 
measured directly and found to vary between 100 /mi and 160 /im. 

The D~ (A+) vertex is fitted in the r-<p plane using all the candidate tracks. The Bg (A^) decay 
vertex is formed by extrapolating the candidate D~ (A+) momentum vector from its decay vertex 
to the intersection with the lepton track. The D~ (A+) decay length is the distance between these 
two decay vertices. The Bg (AP) decay length is found by a fit between the beam spot and the 
reconstructed Bg (A^) decay vertex using the direction of the candidate D~£ + (A+^~) momentum 
vector as a constraint. The two-dimensional projections of the Bg (A^) and Dj (A+) decay lengths 
are converted into three dimensions using the polar angles that are reconstructed from the momenta 
of the B~£ + (A+^~) and D~ (A+). Typical decay lengths for the D~^+ (A+£~) vertex are about 
0.3 cm and the corresponding decay length errors range from about 0.03 cm for the K + K - 7r~, §£~v 
and pK~7r + modes, to about twice this level for the modes which include a A or Kg. 

Additional criteria are used to select D~^ + and A^£~ candidates suitable for precise decay length 
measurements. In channels in which a charm state is fully reconstructed the x 2 °f the charm vertex 
fit is required to be less than 10 (for 1 degree of freedom). Finally, the decay length error of the 
reconstructed Bg (A{J) candidate must be less than 0.2 cm. 

4.4 Results of D~£ + and A+£~ selections 

The invariant mass distributions obtained in each of the D~ decay modes are shown in figure |l[ 
The equivalent distributions for the reconstructed A+ decay modes are shown in figure In each 
case, the fit result overlaid on the histogram is obtained from an unbinned maximum likelihood fit 
to the invariant mass distribution. The results of these fits are summarised in table [l|. 

Each mass fit uses a Gaussian function to describe the signals and a linear parametrization of the 
combinatorial background. In the K + K~7r~ distributions, a second Gaussian is used to parametrize 
contributions from the Cabibbo suppressed decay D — > K + K - 7r~. The mean of this Gaussian is 
fixed to the nominal D~ mass, 1869.3 MeV/c 2 and the width is constrained to be the same as 
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Figure 1: Invariant mass distributions from the four reconstruction channels. In 

each plot, the result of the fit described in the text is overlaid as a solid line. 
The mass ranges used in the decay length fit are shown by the vertical dotted 
lines. 
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Figure 2: Invariant mass distributions from the two A+£~ reconstruction channels. In 
each plot, the result of the fit described in the text is overlaid as a solid line. 
Also indicated are the mass ranges used in the decay length fit. The hatched 
histogram for the [A^ + ]£ _ channel represents the wrong-sign- A combinations: 

[Ai+]r. 



that of the D~ peak. By integrating the tail of the peak due to the D - — > K + K~7r~ decays in 
this D~ signal region, the contamination from this source is found to be negligible. No significant 
peaks are observed in the mass distributions for wrong-sign D~^~ (A^~l + ) combinations in the 
fully reconstructed decay channels K*°K~, 07r~, K~Kg and pK~7r + . The [A£ + ]£~ combinations^] 
include a contribution from A baryons from fragmentation that can be estimated from the wrong- 
sign-A distribution, [A£ + ]£ _ . Studies using simulated data show that the wrong-sign-A distribution 
provides a good representation of these A baryons from fragmentation. The [A£ + ]f~ candidates are 
shown in the plot shaded histogram. 

For each channel, the fitted mass is consistent with the nominal D~ (A+) mass ]l|] and the fitted 
width is consistent with the expected detector resolution. In total, 199 ± 26 D~£ + candidates 
and 145 ± 24 h^£~ candidates are observed. The D~£ + (A^£~) combinations used in the Bg 
(A[J) lifetime fits are selected from regions around the identified invariant mass peaks, including 
a sufficient number of candidates away from the mass peaks to allow an estimate of the lifetime 
characteristics of the combinatorial background. There are 509 for the B^ lifetime fit and 632 for 
the A£ lifetime fit. 



3 The bracketed particles are those that are assigned to be the decay products of a A+. The invariant 
mass requirements on A£ combinations described in sections 4.2 and 4.3, result in a unique assignment of 
the two leptons. 
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Decay 


Signal 


Comb. 


Width 


Fit Range 


Cands. 


channel 


candidates 


fraction 


(MeV) 


(MeV) 


for fit 


K*K 


101±21 


0.46±0.06 


28±5 


-50 - 200 


280 


(f)TT 


53±11 


0.24±0.05 


17±4 


-50 - 200 


114 




8±5 


0.39±0.19 


22±14 


-50 - 200 


21 


<p£ 


37±10 


0.23±0.07 


4 (fixed) 


-8 - 60 


94 


D s 7£+ total 


199 ± 26 








509 


pK7T 


108±22 


0.56±0.06 


21±5 


-50 - 200 


522 


A£ 


37±9 


0.11±0.05 


3.4 ±0.5 


-8 - 80 


69(41) 


k+r total 


145 ± 24 








632 



Table 1: Results of the mass fits to all the signal channels. The second column shows 
the number of candidates in the signal peak. The estimated fraction of combi- 
natorial background is given in the third. The fourth column gives the signal 
peak width found by the fit or the constraint that was used in the fit to the 
invariant mass distribution. The mass range used in the lifetime fit is given in 
the fifth column and the last column gives the number of candidates in this 
mass range. For the [A£ + ]£~ channel the number of wrong-sign-A candidates, 
[A£ + ]£~, is included in brackets. 



4.5 Backgrounds to the — > D g £ + and A£ — ► A+£ signal 

Potential sources of backgrounds to the Bg (A°) signal considered here include decays of other 
b hadrons that can yield a D,r£ + (A+^~) final state or other final states that are misidentified 
as a D~ (A+) hadron. Other sources are D~ (A+) hadrons combined with a hadron that has 
been misidentified as a lepton, and random associations of D~ (A+) hadrons with genuine leptons. 
Finally, there is purely combinatorial background. The various physics backgrounds, and the 
calculation of their contributions relative to that of the signal, are discussed below. 

4.5.1 Physics backgrounds to B° — > D~£ + 

The signal event samples include properly reconstructed D~£ + combinations that do not arise from 
Bg decay. Two decay modes of B° and B + mesons are considered: 

(a) B — > D~DX, D — > £ + v~K (where D is any non-strange charm meson), and 

(b) B — > D~K£ + z/X, where K is any type of kaon. 

For the signal production and decay sequence, the production rate times branching fraction /(b — > 
Bg) • Br(Bg — > D~£ + i/X) = 0.85 ± 0.23% is used Q. Monte Carlo simulations are used to determine 
the selection efficiencies for background modes relative to that of the signal mode. 

For the background, the probability for a bottom quark to form either a B + or B° meson is 
0.378 ± 0.022 |0| each. For the B — > D^K^+z^X mode, the measured branching ratio is less than 
0.009 at the 90% confidence level [GJ. Half of this limit is used as a central value in estimating the 
contribution of this channel, and the range from zero to 0.009 is taken as the uncertainty. For the 
other background mode, it is noted that Br(B — > D S D) = 0.049 ± 0.011 [|l|, which is then corrected 



S 



using Monte Carlo simulation to include the additional contribution from B — » D S DX decays. The 
possibility that all the B — > D S X modes include an additional charm meson is considered as a 
systematic error. The effect on the reconstruction efficiency of D mesons arising from orbitally 
excited D mesons is also taken into account. The total contribution from these two sources of 
backgrounds is estimated to be 11 ± 4%. 



4.5.2 Physics backgrounds to A£ — > A+€ 

The events in the A+ peak may include A.^£~ combinations that do not arise from Ag decay. The 
decay modes considered are B u d — ► A+H C X, H c — » X£~i) and B Uj d, s — ► A,t"X£ _ z/. An estimate of 
the Ag signal contamination from other b baryons is also given. 

As before, the reconstruction efficiencies for these background modes are calculated relative to the 
signal mode from simulated event samples. For the production rate times branching fraction of the 
signal modes, /(b -> Ag) • Br(Ag -» A+£~vX) = 1.35 ± 0.26% § is used. 

To estimate the background from the internal-W decay B — ► A+H C X, H c — > X£~v, the measured 
inclusive branching ratio Br(B — > charmed-baryon X) = 6.4 ± 1.1% jl|] is combined with the 
measurement |2§ Br(B -» A+X)/Br(B -> A~X) = 0.19 ± 0.13 ± 0.04, where B refers only to B 
mesons containing a b quark. It is also assumed that when a A+ is produced, a H c is always 
produced. The average semileptonic branching ratio of the charged and neutral H c (assuming they 
are produced at equal rates) is estimated to be 25 ± 10%. This was obtained from the semileptonic 
branching ratio of the A+, using the theoretical prediction of and the measured lifetimes of 
these baryons After accounting for the relative efficiency, this mode comprises 2.0 ± 1.5% of 
the signal. 

The contribution from the external-W decay, B — > A^X£~u, was estimated using the 90% confidence 
level limit Br(B uc i — » p£~i>X) < 0.16% [^7]]. It is conservatively assumed that this decay always 
proceeds through a A+ and that a A+ is equally likely to decay to a proton or a neutron. Thus, 
Br(B — > A.^£~9X) < 0.32% at the 90% confidence level. The analogous decays of Bg mesons are 
also taken into account by assuming that their branching fraction to A^£~PX is the same as for 
B U)( j and using the hadronization fractions from El . Accounting for the relative detection efficiency 
yields an upper limit for this mode corresponding to 5.7% of the signal. Half of this is taken as the 
central value for this background and the entire range from zero to 5.7% is considered in estimating 
the systematic error. 

Potential contamination of the A.^£~ signal by decays of b baryons other than Ag is also investigated. 
The principle sources of this contamination are from decays of Hb and £]-,. There is some evidence 



for the Hb [28|, which is expected to decay weakly [29]. Theoretical predictions [29] for the Sb mass 
suggest that it is large enough to allow strong decay to a Ag. Accepting this, any non-Ag in the 
signal comes from Hb decays. 

Semileptonic decays of Hb baryons to excited charm-strange baryons that decay subsequently to 
A+ , or non-resonant decays such as Hb — > A+X£ - ^ can contribute to the A.^£~ sample. The level of 
these decays is estimated using the B meson system as a guide. The branching ratio for non-strange 
B decays to X£ + z^, where X is not simply a D or D* meson, is found to be 3.7 ± 1.3%. This value 
is obtained by subtracting Br(B° — > D~ £ + v) and Br(B° — ► D*~£ + v) from Br(B° — > X£ + v) using 
the values from reference ]l]]. The same rate is assumed for the analogous Sb decays mentioned 
above and the conservative assumption is made that in these decays a A+ is always produced. 
It is further assumed that 20% of the weakly decaying b baryons are 3b, based on the relative 
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rates of production of the corresponding light-flavoured baryons in Z° decays Using the above 
branching ratio estimates and reconstruction efficiencies for the signal mode and the modes under 
consideration here, the contribution to the A^£~ signal from the decays of the Hb is estimated to 
be about 1% and is, therefore, neglected. 

The A£ + £~ sample can have additional contributions from the decay of Sb — > H C £~^X followed 
by S c — > Z£ + vX with 3 — > Air. Assumptions about Br(A+ — > AX£+^)/Br(H c — > HX£ + i/) and 
Br(H c — > HX)/Br(A+ — > AX) are required to estimate this background if measurements of A£ |9|, |To| 
and H£ [^8| production are to be used. In the case of the former ratio, this cannot be trivially related 
to the ratio of lifetimes because of Pauli Interference effects between the two strange quarks that 



result from the decay of a H c , but which are not present in the analogous A+ decay |26|]. The latter 
ratio is even harder to predict theoretically, even to the extent of whether one would expect it to 
be greater than or less than unity |}(]]. In estimating the systematic error due to this source, the 
fraction of A£ + £~ candidates due to Sb decays is varied from 10% to 50%. The lifetime will also 
be quoted with a functional dependence on the level of Sb contamination in the A£ + £~ sample, 
fs h /Ai+i~i using the average Hb lifetime of 1. 39^28 P s as measured from E£ correlations 



4.5.3 Other backgrounds 

The D~ (A+) candidate may be a misidentified charm hadron if one or more tracks are assigned the 
wrong particle type. This constitutes an additional source of background which is studied using 
simulated events. For this background, it is found that the invariant mass distribution around 
the D~ (A+) mass is similar to that of the combinatorial background. Such events are therefore 
considered to be included in the combinatorial background fraction. 

The level of background from genuine D~ (A+) particles which are combined with a hadron that 
is misidentified as a lepton can be estimated by fitting the invariant mass spectrum of wrong-sign 
combinations in which the charm candidate and the lepton candidate have the same charge. This 
assumes that random combinations are equally likely to have right and wrong charge correlations. 
For each channel where the charm hadron is fully reconstructed, the wrong sign signal is consistent 
with zero. This is in agreement with what has been found in a related analysis that has greater 



statistical significance [31]. This background source is therefore neglected. 



The background from random associations of a D s (A+) with genuine leptons is estimated using 
simulated data. The contribution is less than one event to each of our samples and is neglected. 

In the two modes in which the charm hadron is partially reconstructed from a semileptonic decay 
channel, D~ — ► <^£ + i/X and A+ — > A^~^X, there are additional backgrounds to consider. These 
include the accidental combination of a ^ (A) , generally produced in fragmentation, with two leptons 
which arise from a semileptonic bottom hadron decay, followed by a semileptonic charm hadron 
decay. In the case of the A^£~ channel, this background is estimated using the observed peak in 
the A invariant mass spectrum for the wrong sign combination formed by an £~ from the bottom 
hadron decay, an £ + from the charm hadron decay and a A (whereas signal would be a A) . Fitting 
the A invariant mass distribution for these wrong sign, [A£ + ]£~ , candidates, yields a contribution 
from this source of 9±5 events. The contribution from random associations of a particle that is not 
from A+ decay with a A and a lepton from A^ decay is found to be negligible. For the analogous 
D~ decay into cp£~VK, there is no wrong sign distribution available, and hence simulated events are 
used to estimate this contribution to be 2.5 ±0.5 candidates. The potential contribution of leptons 
from J/tjj decays, which are then combined with a (j) (either from fragmentation or from a b hadron 
decay) is also estimated using simulated events to be 0.5 ±0.3 candidates. Finally, the contribution 
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from hadrons misidentified as leptons is estimated by selecting events in which the two leptons in 
these modes have the same sign and is found to contribute 2 ± 2 candidates. 

The non-combinatorial background sources mentioned above are expected to contribute a total of 
27 ± 11 events to the D~^ + signal and 16 ± 7 events to the A+£~ signal. The background subtracted 
number of ~D~£ + signal candidates is therefore 

N(B° S -» D~£ + uX) = 172 ± 28 . 
The background subtracted number of A^£~ signal candidates is 

AT(Ag -» A+risX) = 129 ± 25 , 
where no correction has been made for possible Hb contamination. 



5 The B° and A£ Lifetime Fit 

To extract the Bj? (A^) lifetimes from the measured decay lengths, an unbinned maximum likelihood 
fit is performed using a likelihood function that accounts for both the signal and background 
components of the sample. This fit is largely the same as has been used previously in similar 
OPAL measurements || ||. 



5.1 Boost Determination 

For the component of the likelihood function describing the Bg (AP) signals, the Bg (AP) lifetime 
must be related to the observed decay lengths. Since neither channel is fully reconstructed, because 
at least the neutrino produced in the b hadron decay is not reconstructed, it is necessary to estimate 
the b hadron momentum, pb . The probability distribution of a given candidate having a particular 
Bg (Afc 1 ) momentum, B, is estimated on an event-by-event basis in one of two ways, depending on 
the decay channel. 

For the semileptonic D~ decay mode and both A+ decay modes, the technique employed in refer- 
ence H is used. This relies on information from Monte Carlo simulation to estimate the probability 
distribution of b hadron energy, given the observed momentum, p l D , and invariant mass, of all 
the observed tracks in the candidate (i.e., K + K~£~£ + , pK~ir + £~~ or pir~£ + £~). Using a conversion 
factor, R = p 1 d /pb, the relationship of the decay time, t, to the decay length L, can be expressed 
as t = L ■ R • (ms/pl)), where tub is the Bg (A°) mass. The distribution of R was determined using 
simulated data for the signal A+ decay modes and the semileptonic D~ decay mode produced with 
the JETSET 7.4 Monte Carlo 13], using the fragmentation function of Peterson et al. [14]. For 



each decay mode, twelve R distributions were produced covering different ranges of the momentum 
and invariant mass of the D~^ + (A+^~) combination. These distributions are used in the lifetime 
fit to describe the probability, B(pb | Pem^d) a cancn date with a measured p l D and m^, will 
have a particular Bg (A^) momentum. 

For the other D~ decay modes in which the D~ is fully reconstructed, an analytic approach is 



used [||, 31 ] . This approach is not applicable in the case of the A^ decays because of the possibility 
of large polarization effects. The same momentum and invariant mass observables (p l D and m^) as 
the other method are employed, but an exact calculation is used, based on the kinematics of the 
Bg — > D~£ + vX decay rather than Monte Carlo simulation, to calculate the probability distribution 
B(pb I Pj3,m|j). This approach is described in detail in reference |^1J| . 
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5.2 Likelihood Functional Form 



The likelihood function for observing a particular decay length of a Bg (A?) hadron may now be 
parametrized in terms of the measurement error of the decay length, the D~£ + (A+£~) invariant 
mass and momentum, and the assumed lifetime. The functional form of the likelihood is given 
by the convolution of three terms: an exponential whose mean is the Bg (Ag) lifetime, the boost 
distribution obtained from the values of the observed D^rf^ (A+^~) mass and momentum, and a 
Gaussian resolution function with width equal to the measured decay length error. This can be 
expressed as: 



..max 



POO fPj> 

Cf(V I r B ,ai,ph,mh) = dl dp B Q{U \ l,a l L ) B(p B | rfj.mb) V(l \ r B ,p B ) , (1) 

Jo Jo 

where pg ax is the maximum possible energy that the b hadron can have. The function Q is 
a Gaussian function that describes the probability to observe a decay length, L l , given a true 
decay length I and the measurement uncertainty o l L . B is the probability of a particular Bg (Ag) 
momentum for an observed momentum, p l D and invariant mass, m\y of all tracks comprising the 
candidate. V is the probability for a given Bg (A°) to decay at a distance I from the e + e~ interaction 
point. This function is given by: 



r{l | T B ,p B ) = exp 

TbPb 



-I • m B 



(2) 



T B p B 

where T B p B /m B is the mean decay length for a given momentum, p B , mean lifetime, t b , and mass, 
m B , of the B° (Ag). 

As discussed previously, non-combinatorial (physics) backgrounds result from the decay of other 
b hadrons. The likelihood function describing these sources of background is therefore taken to 
have the same form as the Bg (A9) signal, except that the b hadron lifetimes contributing to these 
background samples are fixed to the exclusive world average values Q, weighted appropriately. 
The level of the contributions to this background are set to fixed fractions of the signal, as deter- 
mined in the previous section. The effects of the uncertainty in these fractions on the lifetime are 
addressed as a systematic error. The likelihood function, cf^ M \ describing all sources of D~£ + 
(A+^~) combinations — the B^ (Ag) signal as well as these physics backgrounds — is just a linear 
combination of Cf and the physics background contributions. For the semileptonic channels, the 
lifetime distribution of the backgrounds which include a real <p (A) not from a D~ (A+ ) , is estimated 
from the sideband, in the case of D~, or the sideband and the wrong-sign distribution for the A+ 
mode. For the purposes of determining the lifetime properties of the background, these background 
sources are treated as combinatorial. 

The fit must also account for the combinatorial background present in the D~£ + (A^ - ) sample. 
The functional form used to parametrize this source of background is composed of a positive and a 
negative exponential, each convolved with the same boost function and Gaussian resolution function 
as the signal. This can be expressed as, 

POO „p max 

^ dl J o B dp B G(V I l,*l) B{p B | p^mh.) V bg (L | r+,r 6 - f+,p B ) , (3) 



where 



Pbg(! I rf,T b / 6 + p B ) = fbg ~T~~ ex P 

T bg PB 



I ■ m B 



, n f+\ mB 

+ (1 - / '» ) ^ exp 



(-/) ■ m B 

KglPB 



(4) 



12 



The fraction of background with positive lifetime, f^ g , as well as the characteristic positive and 
negative lifetimes of the background, and T fo ~, are free parameters in the fit. This double- 
exponential shape is motivated by considerations of event topologies that can lead to apparent 
negative decay lengths, even before resolution effects are considered. The background parameters 
are fitted separately for the hadronic and semileptonic D~ (A+) modes, since the lifetime properties 
of the real <p (A) backgrounds may well be different from the purely combinatorial background in 
the other decay modes. 

The background in the event sample is taken into account by simultaneously fitting for the signal and 
background contributions. The probability that a candidate arises from combinatorial background, 
/f omb , is determined as a function of the observed invariant mass of this candidate from the fits to 
the invariant mass spectra shown in figures | andg. 

Thus, the full likelihood for candidate i is: 

C t {V | TB,ai,ph,mh) = (1 - /r mb ) • cf £{M) + /f omb • £? omb • (5) 

In total, four parameters are free in the fit: the Bg (A£) lifetime, and the parameters describing 
the combinatorial background (f^, r^t and r b ~). 



5.3 Lifetime Fit Results 



The fit to the decay lengths of the 509 D~£ + combinations yields T(Bg) = 1.50l ' 15 ps, where the 
error is statistical only. The fit to the 632 A+£~ candidates yields t(A^) = 1. 26^20 P s -Q ^he 
results of these fits are shown in figures [|| and ||. In each figure, the candidates are divided into two 
categories — signal region and sideband region — in order to show the behaviour of the fit when 
the candidate sample consists mostly of signal or mostly of combinatorial background. The signal 
region is defined as the mass region within two standard deviations of the fitted invariant mass (see 
figures [l] and ||). The curves in these figures represent the sums of the decay length probability 
distributions for each event. The figures indicate that the fitted functional forms provide a good 
description of the data for both signal and background. For the Bg (A°) fit, a total \ 2 of 5.2 (9.6) 
is found for the sum of the signal and sideband decay length distributions for 12 (11) bins that 
contain at least five candidates. As was stated earlier, the fits are to unbinned data. 



6 Checks of the Method 



Tests are performed on several samples of simulated events to check for biases in the selection 
and fitting procedures. The first tests involve a simple Monte Carlo program which generates 
decay length data for the signal D~^ + (A^£~) decays and combinatorial background. For each 
signal candidate from a Bg (A^) decay, this simulation generates a Bg (A°) decay time from an 
exponential distribution with the mean set to a known value. The Bg (A[J) momenta are chosen from 
a spectrum based on the full Monte Carlo simulation. The Bg (A[J) decay length is then calculated 
and combined with the momentum to give the true candidate decay time. This is then degraded 
by a resolution function. Physics backgrounds are generated through a similar procedure. Many 

In comparing these two measurements, note that the fractional errors do not scale directly with the 
number of events in the signal because of the larger combinatorial background in the pK~Tr + £~ mode, which 
is the statistically dominant mode for the A° lifetime measurement. 
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Figure 3: Top: The decay length distribution of D,r£ + combinations with an invariant 
mass within the signal region. The unhatched area represents the contribution 
from combinatorial background, the hatched area is the contribution from 
sources of non-combinatorial background and the double-hatched region is due 
to signal from decays of a Bg. Bottom: The similar decay length distribution 
for candidates with an invariant mass in the sideband region. The curves are 
the results of the fit described in the text. 
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Figure 4: Top: The decay length distribution of A+£~ combinations with an invariant 
mass within the signal region. The unhatched area represents the contribution 
from combinatorial background, the (very small) hatched region represents 
non-combinatorial background, not including any Hb contribution, and the 
double-hatched area is due to signal from decays of a A°. Bottom: The similar 
decay length distribution for candidates with an invariant mass in the sideband 
region. The curves are the results of the fit described in the text. 
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fits are conducted over wide ranges of Bg (A°) lifetimes with different levels and parametrizations 
of the backgrounds. The results of these studies show no biases in the fitted Bg (A£) lifetime to 
a level of less than 0.5% and that the statistical precision of the fit to data is consistent with the 
sample size and composition. 

To verify that the D,T^ + (A+£~) selection does not bias the reconstructed sample, lifetime mea- 
surements are made using simulated event samples in which large numbers of the decays of interest 
are produced. In these tests it is found that the mean lifetime of the selected sample of candidates 
is consistent with the lifetime used to generate the sample, indicating no bias in the selection pro- 
cedure. Applying the lifetime fit to the selected samples similarly shows no evidence for a bias to 
within the statistical precision allowed by these samples. This precision ranged from 2% for the 
K + K~7r~ and pK~7r + channels to 4% for the other decay modes. Similarly, applying the selection 
and fitting procedure to a Monte Carlo simulation sample of 4 million hadronic Z° decays, the fitted 
Bg (A^) lifetimes agree with the generated lifetimes to within the statistical power of the sample. 

The lifetime fit is also applied to each decay channel individually. The resulting lifetimes of these 
separate fits are consistent with each other and with the fit to the entire sample. 



7 Evaluation of Systematic Errors 

The sources of systematic error considered are those due to the level, parametrization and source 
of the background, the boost estimation method, possible polarization of b baryons, the beam spot 
determination and possible tracking errors. These systematic errors are summarized in table [| 

The uncertainty in the level of combinatorial background has two sources: the uncertainty due to 
the mass fit to the candidate invariant mass spectra and the statistical variation of the background 
under the invariant mass peak. Background fractions due to one standard deviation variations 
for these two cases are determined and used in the lifetime fit. This produces variations in the 
Bg and A^ lifetimes of ±0.03 ps. The width of the sideband region of the mass spectra, from 
which candidates are selected for use in the lifetime fit, is also varied yielding contributions to 
the systematic errors of ±0.01 ps for the Bg and ±0.03 ps for the A^ lifetimes. Using a quadratic 
function to describe the mass distribution of the combinatorial background has a negligible effect 
on the fitted lifetimes. 

The effect of the uncertainty in the level of the non-combinatorial backgrounds to the Bg and A^ 



signal is estimated by varying these background levels over the ranges described in section 4.5. This 
produces a ±0.01 ps variation in the Bg lifetime. The systematic shift due to a Hb contamination 
of 30 ± 20% in the A° sample is ±0.02 ± 0.02 ps. The shift in the central value of the measured 
lifetime as a function of Hb contamination is ±0.08 • fs h /Ai+£- P s - All other non-combinatorial 
backgrounds mentioned in section 4^5 contribute an additional error of ±0.01 ps to the A^ lifetime. 



The b hadron lifetimes used for these backgrounds are also varied within their measured errors |lj. 
The resulting change in the fitted lifetimes is less than ±0.01 ps. 

The total systematic error associated with the description of the combinatorial and physics back- 
grounds is therefore ±0.03 ps (±0.05 ps) for the measured Bg (A^) lifetimes. 

The effects of uncertainty in the b hadron fragmentation are estimated slightly differently for the 
two boost estimation methods employed. For the hadronic D~ decay modes, the estimated Bg 
energy spectrum used by the boost estimation procedure is varied within the measured limits of 
the average b hadron energy Similarly, in generating the Monte Carlo events used to estimate 
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the b hadron boost for the two A+ decay modes and the (f>£~ mode, the average b hadron energy was 
varied by the same amounts as used above. These variations yield changes in the fitted lifetimes of 
±0.02 ps. The effect on the lifetime of varying the mass of the by ±50MeV/c 2 about its central 
value of 5641 MeV/c 2 |J is ±0.01 ps. The effect of a 2 MeV/c 2 uncertainty in the mass of the B° |l| 
results in a change of less than 0.01 ps in the measured Bg lifetime. 

In the Monte Carlo events used for the boost estimate, the A^ was assumed to be unpolarized. 
However, in the Standard Model, b baryons can retain up to the full longitudinal polarization of 
—0.94 from the b quark. A variation from to —0.94 polarization produces a change of ±0.06 ps. A 
recent measurement of the polarization |3^] is used to correct the lifetime extracted using the decay 
length fit. This yields a correction of ±0.014]to'oi4 ps, where this error results from the precision of 
the polarization measurement. The effect of the choice of form factor used to describe the energy 
transfer from the A^ to the A+ has also been investigated. The use of the alternative form factors 
of reference |l7| produces a negligible change in the fitted lifetime. 

The average interaction point of the LEP beams in OPAL is used as the estimate of the production 
vertex of the Bg and A^ candidates. The mean coordinates of the beam spot are known to better 
than 25 fim in the x direction and 10 /um in y. The effective r.m.s. spread of the beam is known 
to a precision of better than 10 /j,m in both directions. To test the sensitivity of r(Bg) and t(A^) 
to the assumed position and size of the beam spot, the coordinates of the beam spot are shifted 
by ±25 ^m, and the spreads are changed by ±10 /im. The largest observed variation in r(Bg) and 
r(A[J) is 0.01 ps which is assigned as a systematic error to both measurements. 

The effects of alignment and calibration uncertainties on the result are not studied directly but are 
estimated from a detailed study of 3-prong r decays [p4fl , in which the uncertainty in the decay 
length due to these effects is found to be less than 1.8% for the data taken during 1990 and 1991 
and less than 0.4% for later data. This corresponds to an uncertainty on r(Bg) and t(A^) of 
0.01 ps. The potential for incorrect estimation of the decay length error is addressed by allowing 
an additional parameter in the lifetime fit which is a scale factor on the estimated decay length 
error. This parameter is found to be consistent with unity. This procedure changes the Bg lifetime 
by less than 0.01 ps and the A° lifetime by — 0.02ps. 



Source 


r(B s ) ) correction (ps) 


t(A^) correction (ps) 


background (excl. E^) 


0.00 ±0.03 


0.00 ±0.05 


Hb background 




±0.02 ±0.02 


uncertainty in boost 


0.00 ±0.02 


0.00 ±0.02 


polarization 




±0.01 ±0.02 


beam spot 


0.00 ±0.01 


0.00 ±0.01 


alignment errors 


0.00 ±0.01 


0.00 ±0.01 


total 


±0.04 


±0.03 ±0.06 



Table 2: Summary of systematic corrections and uncertainties on the Bg and A^ life- 
times. 



8 Conclusion 

The decay channels B° -> D-£ + vX and Ag -> JS+£~vX, where D~ decays to K*°K", 0vr~, K"Kg 
or (j)£~vX and A+ decays to pK~7r + or A£ + vX have been reconstructed. From almost 4.4 million 
hadronic Z° events recorded by OPAL from 1990 to 1995, a total of 172 ± 28 such candidates are 
attributed to Bg decays and 129 ± 25 such candidates are attributed to A^ decays. 
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The Bg lifetime is found to be 



r(B s °) = 1.50l°il±0.04ps, 

where the first error is statistical and the second systematic. As predicted by theoretical calcula- 
tions [||, this result is consistent with the observed value for the B° lifetime Q. This is also in 
agreement with other measurements of the Bg lifetime || fj|. The above value of r(Bg) has been 
combined with the OPAL measurement of the Bg lifetime in which only a D~ candidate is recon- 
structed H which yields r(Bg) = l^+Q^stat^Q^syst) ps. Taking the correlated statistical and 
systematic errors into account, the average of these two measurements is found to be 1.57±0.14ps. 

The measured A^ lifetime, is 

r(A°) = (1.27 + 0.08 • f Sb/M+e -)±ail ± 0-06 ps , 

where the dependence of the fitted A^ lifetime is given in terms of the fraction, f~, h /Ae+£- j of the 
A£ + £~ candidates that are due to Eb decays, the first error is statistical and the second systematic. 
Assuming Hb contamination of 30 ± 20%, the A^ lifetime is 

r(A°) = 1.29±°;i±0.06ps. 

The lifetime using the A£ + £~ sample alone is found to be O.SSIq'^Iq'^ ps, when the Hb content 
is taken to be 30 ± 20%. For the pK~ir + £~~ sample alone, the Eb content is estimated to be only 
about 1% and the lifetime of this sample is 1. 361^23 -° 0-06 ps. These are consistent with other 
recent measurements [p], 10] which have tended to be lower than the B° meson lifetime in qualitative 



agreement with the predictions of || 
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